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Abstract 

Imposing a constraint of the sound-speed profile determined from helioseismology and updating the 
microphysics, we have revised our seismic solar model, constructed with the assumption of a homogeneous 
metal abundance distribution, and have shown that the theoretically expected neutrino fluxes are still 
significantly more than the observations. With the same sound-speed profile constraint, we also constructed 
solar models with low metal abundance in the core, and evaluated the neutrino fluxes of these models to 
see if nonstandard solar models with a low metal core can solve the solar neutrino problem. Some of 
these models are in agreement with the Homestake data, the Super-Kamiokandc data, and the sound- 
speed profile simultaneously, but none of these satisfy both the neutrino flux data, including GALLEX and 
SAGE, and the helioseismically determined density profile. 

Key words: neutrinos — Sun: abundances — Sun: helioseismology — Sun: interior — Sun: oscilla- 
tions 



1. Introduction 

Various nonstandard solar models have been proposed 
to explain the deficit of neutrino fluxes from the Sun com- 
pared with the theoretical expectation based on the stan- 
dard evolutionary models of the Sun. Solar models hav- 
ing a core with low heavy clement abundance, Z, arc such 
examples (Joss 1974; Levy, Ruzmaikina 1994; Jcffery et 
al. 1997). It is expected that a low Z in the central re- 
gion of the Sun would lead to low opacity there, result- 
ing in a lower central temperature; such models are ex- 
pected to be in better agreement with the observed neu- 
trino fluxes. However, nonstandard solar models have of- 
ten been criticized: either their properties are in substan- 
tial disagreement with seismic data, or their sound-speed 
profile disagrees significantly with the seismically deter- 
mined sound-speed profile of the Sun (e.g. Bahcall et al. 
1997; Gough 1999). For example, comparing the observed 
data and the model predictions of the frequency difference 
between the n-th overtone p-mode of the low-degree I and 
the (n — l)-th overtone p-mode of the degree 1 + 2, which 
is sensitive to the structure of the solar central region, 
Guenther and Demarque (1997) claimed that low-Z cores 
extending beyond O.OBMq should be ruled out because, 
otherwise, the model predictions would be in disagreement 
with the observed data and, furthermore, there would be 
no chance that the low-^ core models could reduce the 
neutrino capture rates below ~4.87 SNU and ~119 SNU 
for the Homestake neutrino detector using ^^Cl and the 
GALLEX and SAGE detectors using ^^"^Ga, respectively, 
while the detected capture rates are ~2.56 ± 0.23 SNU 
and ~72 ± 8 SNU, respectively (A SNU is defined to be 
10""^^ interactions s~^ per target atom). Such criticisms 
are, however, not necessarily fair, because those models 



were constructed without fine tuning to agree with the 
helioseismic data. 

The recent observations of solar oscillations provide us 
with a large number of very accurate eigenfrequencies of 
the Sun. The relative errors in the frequency measurement 
are as small as 10~^. With the inversion of these frequen- 
cies, we now know the sound-speed profile in the Sun, c(r), 
within errors of a few tenths of a percent. Imposing the 
constraint of the thus-obtained sound-speed profile, the 
profiles of temperature, density, pressure, luminosity, and 
chemical composition can be deduced (Shibahashi 1993; 
Shibahashi, Takata 1996). 

It should be stressed here that the sound-speed profile 
has been very precisely determined from helioseismology, 
and that the difference between the profile determined 
in this way and that of any standard evolutionary solar 
model is demonstrably larger than the observational er- 
ror. This means that we can also determine the profiles of 
thermal quantities in the Sun by using the observationally 
based highly accurate profiles for the acoustic quantities 
and adopting the same information of the microphysics 
(equation of state, nuclear reaction rates, and opacity) as 
used in constructing standard evolutionary solar models. 
The benefits of building such a solar model are as follows 
(Shibahashi 1993, 1999; Takata, Shibahashi 1998a, here- 
after TS98a): 

• We can construct a model of the present-day Sun 
without any assumption about the evolutionary his- 
tory of the Sun. For example, we need not worry 
whether the ■^He induced g-mode instability and/or 
mass loss affects the evolution of the Sun. 

• The location of the base of the convection zone, 
''conv, can be estimated from the fact that the sound- 
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speed profile changes rapidly there (Gough 1984). 
Knowing the depth of the convection zone allows 
us to treat the radiative core only without worry- 
ing about the treatment of convection, which is not 
well-described theoretically (TS98a). 
• The model is naturally consistent with the helioseis- 
mically determined sound-speed profile, while the 
evolutionary solar models are not necessarily so. 

The main aim of this paper is to discuss our construc- 
tion of solar models with an assumption of low Z in the 
core with the imposition of a constraint of the sound-speed 
profile, and to evaluate the neutrino fluxes of these models 
to see if nonstandard solar models with a low-Z core can 
solve the solar neutrino problem. Solar models based on 
helioseismology were constructed with an assumption that 
Z is homogeneous throughout the Sun (TS98a). Since 
then, the nuclear cross sections have been systematically 
recompiled and updated. Before discussing low-Z core 
models of the Sun, in section 3 we revise the seismic so- 
lar model with a homogeneous Z by adopting updated 
input physics. Then, in section 4, with the same input 
physics, we construct seismic solar models having a low-Z 
core and evaluate their neutrino fluxes. The preliminary 
results were presented at the SOHO-lO/GONG-2000 joint 
workshop (Watanabe, Shibahashi 2000). We describe here 
our method in more detail, and present the latest results, 
including an estimate of the errors in the model properties 
and in the neutrino fluxes. 

2. Equations of the Seismic Solar Model 

We solve the basic equations governing the radiative 
core of the Sun with the imposition of information from 
helioseismology (TS98a, Takata. Shibahashi 1998b, here- 
after TS98b, 1999; Shibahashi et al. 2000). They are the 
continuity equation, 

dMr 



e = e{P,T,Xi). 



(7) 



dr 



the hydrostatic equation, 

dP _ GMrP 

dr ' 
the energy equation, 

dLr 2 
— — = Airr pe, 
dr 

and the energy transfer equation, 
dT 3 Kp Lr 
dr 



(1) 



(2) 



(3) 



(4) 



4ac T3 47rr2 ' 
where the symbols have their usual meanings. 

In addition to these differential equations, we need aux- 
iliary equations: the equation of state, the equations for 
the opacity and the nuclear reaction rates: 

p = p{P,T,X,), (5) 

n^K(P,T,X,), (6) 

and 



These equations link the thermal quantities and the chem- 
ical abundances. In making standard evolutionary mod- 
els, the chemical abundance distribution is obtained by 
following the time evolution. 



dX, 
dt 



dX, 
dt 



nuclear 



dX, 
dt 



(8) 



diffusion 



where the nuclear reactions and the diffusion processes are 
taken into account. In the standard scenario, the chem- 
ical composition is assumed to be uniform at zero-age; 
dXi{t ~ 0)/dr ~ 0. In many nonstandard evolutionary 
models, some other time evolution processes are intro- 
duced and/or the initial conditions at t = are different 
from the standard case. 

On the other hand, in constructing seismic models, we 
do not follow the time evolution. Instead, we impose the 
seismically determined sound-speed profile as a constraint 
on the model. Note that, if we distinguish only the hy- 
drogen (^H) and helium f'He) separately as X and Y, 
respectively, and treat all other elements collectively as 
heavy elements, Z, then the sound-speed, as one of the 
thermodynamical quantities, is a function of two other 
thermodynamical quantities (P and T) and X and Z, 



c{P,T,X,Z)=CoUr)- 



(9) 



This inversely relates the hydrogen abundance, X, with 
the pressure, temperature, metal abundance, and the 
sound-speed, 

x = x(p,r,z,cobs). (10) 

Since X is given by equation ( |lO| ) with the helioseismi- 
cally determined sound-speed, the density, the opacity, 
and the nuclear reaction rate are, in turn, given in terms 
of (P,T,Z,Cobs) by equations (5)-(7). 

If Z is given, then with the help of equation (10), the 
basic equations (l)-(4) can be solved with the proper 
boundary conditions. In principle, the helioseismically de- 
termined density profile. pobs(''), can be used as another 
constraint in constructing a seismic solar model. This 
constraint together with the sound-speed constraint de- 
termines the heavy elements' abundance profile, Z, at a 
given r as a part of the solution. Such attempts have in- 
deed been carried out (TS98b). They have, however, not 
yet succeeded in obtaining a reasonable Z-profile in this 
way, since the dependence of the equation of state upon 
Z is so weak. We should also note that the formal error 
of Pobs('') is much larger than that of Cobs(''). In this pa- 
per, therefore, we impose only the sound-speed profile as 
a constraint in constructing a seismic solar model, and we 
assume a certain Z-profile a priori. 

In section 3, Z is assumed to be homogeneous for sim- 
plicity, following TS98a, while the microphysics is up- 
dated. In section 4, we assume that in the central core 
the heavy elements' abundance, Z, is lower than the outer 
part: low-Z core models. The schematic Z-profiles are 
shown in figure 1. A typical Z-profile of a standard evo- 
lutionary model with gravitational settling is shown in 
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figure 1 c (e.g., Bahcall et al. 1998, hereafter we refer to it 
as BP98). TS98b treated a rectilinear Z-profile and their 
Z-profile is just like figure 1 d. 

The boundary conditions at the center are trivial: 

Mr = and L, = at r = 0. (11) 

The location of the base of the convection zone is es- 
timated from the inverted sound-speed profile (Gough 
1984). We adopt the seismically determined depth of the 
convection zone, r = rconv, and set the outer boundary 
conditions at the base of the convection zone (TS98a): 



V 



ad 



ad 



WnacGMrT^ 

at T ^conv • 



and Lr ~ Lp 







(12) 

This means that we do not need to care about the con- 
vective heat transport, which has theoretical uncertain- 
ties. Moreover, chemical homogeneity in the convection 
zone requires that Z/X at the base of the convection 
zone, (Z/ X)convi should be identical with the value at 
the photosphere, {Z/X)surf, which is determined spectro- 
scopically. 

Table 1 summarizes the microphysics and parameters 
of the seismic solar models adopted in the present work. 
We adopt the updated microphysics which was also used 
by BP98. Table 2 summarizes the differences in some im- 
portant nuclear cross-section factors between TS98a and 
this work (Adelberger et al. 1998). As for the helioseis- 
mically determined sound-speed profile, we adopt Basu's 
(1998) results obtained by inversion of the MDI veloc- 
ity data of the first 144 days. The formal error of the 
sound-speed inversion is smaller than 0.05% in the region 
0.2 < r/i?0 < 0.8, and does not exceed 0.3% in the entire 
region. In the central core {t/Rq < 0.05), we extrapolate 
the inverted sound-speed profile and its error level to the 
center. 

3. Structure of the Seismic Solar Model with 
Homogeneous Z 

We first construct a seismic solar model by assuming 
that Z is homogeneous. The methodology is exactly the 
same as that adopted by TS98a. The properties of the 
homogeneous Z model are summarized in table 3. 

The effect of various uncertainties in microphysics upon 
the seismic model and the theoretically expected neutrino 
fluxes was investigated by a Monte Carlo simulation. For 
example, as for each of the nuclear cross sections, we con- 
structed one hundred sets of seismic models by super- 
imposing Gaussian noise corresponding to the l-cr level 
uncertainty of each one of the astrophysical S'-factors, 
while keeping other microphysics unchanged so as to iso- 
late the effect of each microphysics. As for the inverted 
sound-speed profile, we also constructed one hundred sets 
of seismic solar models by adding Gaussian noise to the 
most likely sound-speed at every step of At/Rq = 0.01. 
In this process, we interpolated the sound-speed at the 
other mesh points smoothly. We estimated the effect of 
the uncertainties in opacity by constructing an additional 



one hundred sets of seismic models with opacities having 
Gaussian noise with 5% amplitude added in the whole 
range of density, temperature and chemical compositions 
at all mesh points. To see the effect of the equation of 
state, we constructed a seismic solar model with the ideal 
gas law, and compared the model structure with that con- 
structed with the OPAL equation of state. Table 4 sum- 
marizes the effect of the uncertainties in microphysics, in- 
vestigated in this way, upon the theoretically expected 
neutrino fluxes, the central density, pc, and the helium 
abundance at the surface and the base of the convection 
zone, Fconv, of the seismic solar model. 

It is clear from this table that the influences of 
the S'(0)-factors of ^ll{p,e+i^^)^ll-, 3He(4He,7)^Bc-, and 
''Be(p,7)®B-reactions are crucial for the neutrino fluxes. 
The increase in the S'ii(0)-factor makes the nuclear energy 
generation more efficient. Since the luminosity should be 
fixed; this means that the temperature and density near 
the center should become lower. Since the sound-speed 
(c^ = TiP/p « TiRT/fi) is also fixed, the decrease in the 
temperature, Tc, and in the density, pc, leads to an in- 
crease in the hydrogen abundance, Xc, and a decrease 
in the pressure, Pc, respectively. Decreases in pc and 
Tc make the pp-II and the pp-III reactions less efficient. 
Thus, the ^Be- and the '^B-ncutrino fluxes decrease and, 
hence, the neutrino capture rates for both the chlorine 
experiment and Supcr-Kamiokandc become smaller. The 
decrease in the ''Be- and the ^B-neutrino fluxes dominates 
over a slight increase in the pp-neutrino flux; the neutrino 
capture rate for the gallium experiments also decrease. 
The increase in the S'33(0)-factor makes the pp-I reaction 
more efficient, and thus the contributions of the pp-II- and 
pp-III reactions to the energy generation become smaller. 
Hence, the ^Be- and the ^B-neutrino fluxes caused by the 
pp-II and the pp-III reactions, respectively, decrease, and 
then the neutrino capture rates for the chlorine experi- 
ment and the Super-Kamiokande decrease. Decreases in 
the ''Be- and the ^B-neutrino fluxes also lead to a de- 
crease in the neutrino capture rate for the gallium exper- 
iments. On the other hand, an increase in the £'34(0)- 
factor makes the branching ratio of the pp-II and pp-III 
reaction to the pp-I reaction larger. Consequently, the 
^Bc- and the ®B-neutrino fluxes increase, and the neu- 
trino capture rates for all of the neutrino detection exper- 
iments increase. The increase in the 5i7(0)-factor makes 
the branching ratio of the pp-III reaction to the pp-II re- 
action larger. Hence it makes the ^B-neutrino flux and the 
neutrino capture rates for the chlorine experiment and the 
Super-Kamiokande particularly larger. For the neutrino 
capture rate for the gallium experiments, the uncertainty 
in the neutrino cross-section is also crucial because the 
uncertainty is larger than that for the chlorine detector. 
On the other hand, the mass ratio of metal to hydro- 
gen at the surface, {Z/X)surf, the depth of the convection 
zone, Tconv/Re, the opacity, and the equation of state are 
crucial for the surface helium abundance, l^conv Among 
them, {Z/X)sur{ is crucial because it is related to yconv 
directly and its uncertainty is large. 

A comparison of the present result with observations 
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and some other solar models are summarized in table 
5 and depicted in figure 2. Note that errors of TS98a 
and TS98b are rough estimations. The seismic solar 
model presented by TS98a was constructed with the same 
assumption of Z as in the present case. We thus re- 
gard the present model as a revised version of TS98a. 
While the sound-speed profile is replaced with the up- 
dated one, the difference is so small that this replacement 
does not induce a significant difference in the structure of 
the model. Rather, the difference in the model structure 
comes mainly from the updated input microphysics. As 
can be seen in table 5, the neutrino-capture rates actually 
detected in the current experiments are still significantly 
less than predictions based on the present seismic solar 
model. 

Figure 3 shows the relative differences in the sound- 
speed profile and in the density profile between the present 
seismic solar model and a direct inversion of the oscilla- 
tion data (Basu 1998). For a comparison, the differences 
between the standard evolutionary model (BP98) and the 
direct inversion are also shown. Since we imposed a con- 
straint of the inverted sound-speed profile in construct- 
ing the seismic solar model, the relative difference in the 
sound-speed of the seismic model shown in the left panel 
is exactly equal to the amount of the uncertainty in the in- 
verted sound-speed profile, itself. Note that the difference 
between the sound-speed profile of the evolutionary stan- 
dard model (e.g., BP98) and the helioseismically inverted 
profile is larger than the 1-cr level of uncertainty. The seis- 
mic model is consistent with the inverted density profile 
only within the l-a level, while the difference between the 
density profile of the evolutionary standard model (BP98) 
and the inverted profile is larger than the 2-cr level near 
the base of the convective zone. 

The uncertainty in the density profile of the seismic so- 
lar model is, however, larger than the 1-cr level uncertainty 
of the inverted result. Note that the origin of each uncer- 
tainty is different. The uncertainty of the inverted result is 
caused by the measurement errors of solar oscillations. In 
order to see what is the main cause of the large uncertainty 
of the seismic solar model, we carried out a Monte Carlo 
simulation while changing only one input microphysics, 
and keeping all others unchanged to isolate the effect of 
each microphysics. Figure 4 shows the uncertainty in the 
density profile, caused only by each input microphysics. 
The influence of the highly accurate sound-speed profile 
is very small. From this simulation, it is concluded that 
the main causes are the nuclear cross-sections of the pp- 
reaction [S'ii(O)] and of the ^He-|- "'He-reaction [S'34(0)]. 
Since energy generation in the Sun is mainly controlled 
by 511(0), it is most crucial for the solar structure: den- 
sity profile. Although 534(0) is less important than 5ii(0) 
for the solar structure, its uncertainty is much larger. 

Figure 3b shows that the uncertainty in the density pro- 
file of the seismic model is very small around r = 0.2Rq. 
Since the total luminosity is fixed, with the increase of 
5ii(0), which is most important for the density profile, 
the density in the core decreases to compensate for the 
increase in the nuclear reaction rate. On the other hand. 



since the total mass is fixed, the density outside the core 
increases to compensate for the decrease of the density 
in the core. Because these opposite reactions of the den- 
sity against the input physical parameter balance around 
r = O.2i?0, the uncertainty in the density of the seismic 
model is very small around there. 

4. Structure of the Seismic Solar Model with a 
Low-Z Core 

Solar models having a low-Z core are expected to lead to 
low neutrino fluxes. One possible scenario for the forma- 
tion of a \ow-Z core is to assume that dust and heavy ele- 
ments remaining in the disk around the Sun after the for- 
mation of the planets were accreted onto the Sun so that 
heavy elements would be accumulated more in the enve- 
lope than in the core (Joss 1974; Jeffery et al. 1997). In an- 
other unconventional scenario, heavy elements locked up 
in grains are assumed to be segregated from the hydrogen 
and helium gas in the pre-solar nebula (Levy, Ruzmaikina 
1994). In order to see how much the neutrino fluxes can 
be reduced by introducing a low-Z core, we constructed 
seismic solar models having a low-Z core with the impo- 
sition of the sound-speed profile. We assume here that 
the Z-profile is a step function of r; that is, Z(r) is a 
certain constant, Z^, in the region < r < rf, where rf 
is a parameter, while Z is another constant in the region 
ff'£fl£ fconv, so that Z/X matches its surface value, 
{Z/X)sui-i (figure 1 b). The parameter Zc ranges from 
0.0001 to 0.020. We constructed the models so that their 
sound-speed profiles would be consistent with the seismi- 
cally determined profile. Table 6 shows the properties of 
low-Z core models with various values of rt in the case 
of Zc = 0.0001. Note that the model with rf = is the 
same as the homogeneous Z model discussed in the previ- 
ous section. In this table, Mcorc denotes the mass of the 
low-Z core [Afcoro = Mr{r = rf)]. 

The properties of the low-Z core models can be ex- 
plained as follows. With an increase in the low-Z core 
size, the opacity in the core decreases, resulting in a de- 
crease of the central temperature. Since the sound-speed 
profile is fixed by the observation, this requires a decrease 
in the mean molecular weight to compensate for a decrease 
in the temperature, and hence an increase in the hydrogen 
abundance there. The total luminosity is also fixed, which 
requires an increase in the density in the core to compen- 
sate for the decrease in the nuclear reaction rate due to 
the decrease in temperature there [cf. equation (3)]. The 
constraint of the sound-speed leads to an increase in the 
pressure in the core because of the higher density. These 
tendencies are common in the seismic solar models with 
a low-Z core, and become conspicuous with a decrease in 
the value of Zc. 

Figures 5 a-c show the dependence of the neutrino cap- 
ture rates of the seismic models upon the values of Zc and 
Afcorc- For a comparison, the detected neutrino capture 
rates arc also indicated. The neutrino capture rates for 
the chlorine detector (Homestake) and the ^B-ncutrino 
flux (Super-Kamiokande) decrease substantially with an 
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increase in the low-Z core size. 

The effect of various uncertainties in microphysics upon 
the model properties and the thcorcticaUy expected neu- 
trino fluxes was investigated by a Monte Carlo simulation 
in the same way as that described in the previous sec- 
tion. The results are summarized in table 7 for the case 
of Zc = 0.0001; also, the thus-determined error bars for 
the neutrino capture rates, the *B-neutrino flux and the 
central density are depicted in figure 5. 

Let us focus attention on those models with = 0.0001 
(the continuous line). As can be seen in figure 5 a, the 
models with Mcorc > O.IM© do not contradict the result 
of the Homestake experiment if we tolerate the differ- 
ence within the 1-a level. In the same sense, the Super- 
Kamiokande experiment is consistent with those models 
having O.OSMq < Mcore < O.12M0. That is, low-Z core 
models with Mcore ~ O.IMq are favorable for explaining 
the Homestake and the Super-Kamiokande experiments 
simultaneously. Besides that, such models are naturally 
consistent with the seismically determined sound-speed 
profile. These results are not in agreement with the con- 
clusion of Guenther and Demarque (1997, hereafter we re- 
fer to it as GD97) , who investigated the low-Z core model 
of the Sun by following its evolution. There are differ- 
ences between the present approach and GD97's in the 
concept and methodology of making solar models. While 
the models of GD97 are constructed following the stel- 
lar evolution, the present models are not. In this sense 
the present approach has more freedom. We do not try 
to explain how our seismic models can be realized in the 
evolutionary process, and such an investigation is beyond 
our scope. While the present models were constructed 
with fine tuning to agree with the helioseismic data, those 
of GD97 were not. It also should be noted that GD97 
adopted almost the same input physical parameters with 
TS98a. Updating them in GD97's approach may lead to 
lower neutrino fluxes (cf. figure 2 and table 5), and the 
difference between the present result and GD97's would 
become smaller. Our results demonstrate that certain 
low-Z core models of the Sun can be in agreement with 
the Homestake and the Super-Kamiokande data and he- 
lioseismology. Therefore, criticisms based on evolution- 
ary models are not necessarily fair. However, even the 
present seismic solar models are substantially in disagree- 
ment with the neutrino capture rates for the gallium de- 
tectors, GALLEX and SAGE. 

In the present approach, we did not use the seismically 
obtained density profile as a constraint in making seismic 
solar models. We should compare the density profile of 
the models with the seismically inverted profile. In the 
case of Zc = 0.0001, as shown in figure 5 d, the central 
density of the model becomes higher with an increase in 
the core size, and only those models with Mcore < 0.02Mq 
are consistent with the seismically inverted density (Basu 
1998). To see the effect of parameter Zc, we adopted 
various values of Zc and constructed seismic models of the 
Sun with a low-Z core in a similar way. However, no model 
is consistent with all of the neutrino experiments and the 
inverted density profile. From these results, we conclude 
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that none of the low-Z core models are simultaneously 
consistent with all of the neutrino experiments and the 
helioseismically inverted density profile. 

In the case of Zc > Zgurf (~ 0.018), the temperature at 
the central region becomes higher than that of the homo- 
geneous Z model, and hence deviations from the detected 
solar neutrino fluxes become larger, while the density at 
the central region becomes less inconsistent with the seis- 
mically inverted proflle. 

5. Summary and Conclusion 

We have revised the seismic solar model with a con- 
straint of the sound-speed profile, and the latest input 
microphysics and physical parameters, assuming that Z 
is uniform, and also have estimated the uncertainty of 
the seismic model thoroughly by a Monte Carlo simula- 
tion. Though the seismic solar model thus constructed 
is naturally consistent with the helioseismic data, the 
expected neutrino fluxes based on this model are still 
higher than the experimental data. This result implies 
that a modiflcation of the solar model does not seem to 
be able to solve the solar neutrino problem, leaving the 
neutrino-oscillation explanation as an attractive possibil- 
ity. The uncertainty in the seismically determined sound- 
speed profile is so small that it does not induce substan- 
tial uncertainties in the physical properties of the model. 
The uncertainties in the nuclear cross-section [especially 
5*11(0), 534(0), and 5*17(0)] have more crucial influences. 

We have also constructed seismic models with an as- 
sumption of a low-Z core, and evaluated the neutrino 
fluxes of these models to see if these nonstandard solar 
models can solve the solar neutrino problem. Although it 
turns out to be possible to explain the neutrino fluxes for 
both the Homestake and Super-Kamiokande experiments 
and the sound-speed proflle simultaneously by the low-Z 
core extending ^ O.OlAfcorc/M0, it is still impossible to 
explain all of the data, including the neutrino flux mea- 
surements based on the gallium detectors, by this model. 
The density in the core becomes higher with a decrease of 
Zc, and the model becomes less consistent with a seismi- 
cally inverted density profile. 

Therefore, we conclude that the low-Z core models of 
the Sun cannot explain all of the solar neutrino detection 
experiments and the helioseismically inverted sound-speed 
and density profiles. 
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Table 1. The Microphysics and physical parameters of the seismic solar models. 



Physics /Parameter 




Reference 


Equation of state 


OPAL 


Rogers et al. 1996 


Opacity 


OPAL 


Iglesias, Rogers 1996 


Nuclear cross-section 




Adelberger et al. 1998 


Screening effect 


Weak screening 


Salpeter 1954; Gruzinov, BahcaU 1998 


Neutrino cross-section 




Bahcall 1997; BahcaU et al. 1996, 1998 
Garcia et al. 1991 


{Z/X)su-cf 


0.0245 


Grcvesse, Noels 1993 


Luminosity Lq 


3.844(1 ±0.004) X lO^^erg s^^ 


Bahcall, Pinsonneault 1995 


Sound-speed profile 


SOHO data 


Basu 1998 


^conv 


(0.713 ±O.OOl)i?0 


Basu, Antia 1997 



Table 2. Comparison of some important nuclear cross-section factors between TS98a and this work. 





TS98a 


This work 




Reference 


BahcaU, Pinsonneault 1995 


Adelberger et al. 


1998 




S'(O) (keV barns) S' (0) (barns) 


S'(O) (keV barns) 


S'{Q) (barns) 


iH(p,e+:/e)^H .. 
3He(3He,2p)4He 
3He(''He,7)^Be . 
^Be(p,7)8B .... 
i'*N(p,7)i50 ... 


3.89(l±0.011)xl0-22 4.52x10-24 
4.99(l±0.06)xl03 -0.9 
0.524(1±0.032) -3.1 x 10"'' 
0.0224(1±0.093) -3 x 10-^ 
3.29(1±0.12) -5.91x10-3 


4.00(l±0.007l^'^^^) X 10-22 
(5.4±0.4)xl03 
0.53±0.05 
0.019l°:g°^(la)tg:r4(3a) 
3.5t?:^(la)t^:^(3a) 


4.48x10-24 

-4.1 
-3.0 X lO-'' 
-1.3 X 10"^ 
-1.28 X 10-2 



Table 3. Properties of the seismic solar model with homogeneous Z. 



Quantities 

Tc(10^ K) 

Pc (10" dyn cm-2 ) 

Pc (g cm-3 ) 

X, 



M^ouv/Mq 

Teonv (106 K) 

Pconv (10^3 dyn cm-2 ) 

Pconv (g Cm-3 ) 

A^conv 

Y 

^ conv 

^conv 

pp V flux (10^° cm-2 S"^ ) . 
pep v flux (10^ cm-2 g-i -j ^ 
hep V flux (10'^ cm-2 g-i ^ ^ 
''Be V flux (10^ cm-2 g-i ^ ^ 
^B ly flux (10^ cm-2 g-i ) ^ ^ 
i^N 1/ flux (10^ cm-2 g-i ) ^ 
1^0 flux (10*^ cm-2 g-i j ^ 
i^F V flux (10^ cm-2 g-i ) 
v capture rate for CI (SNU) 
ly capture rate for Ga (SNU) 



Most likely values 

HO. 005 



1.561 _o.oo9 

9 Q7Q +0.031 
z.0(0 _o.o49 

156.0 ^2 ° 

n •^•ifi'i 4-0.0058 
U.OOOO _0.0035 

6437 +0 0035 

U-D^-J' -0.0058 



-Q-Q264-FinroDr- 

V.U/.U'i _0.0013 

2.19 ±0.01 

r 71 -1-0.07 

"J-'^ -0.06 

0.190 ±0.001 

n TiP.^ +0.0036 
V.I 000 _o.o037 

2455 +0 0042 

U.Z400 _o.0041 

0.0180 ±0.0004 



5.98 ±tt 
1.44 ±0.02 

2.11 ±0.06 

+0.39 
-0.43 

+ 1.04 
-0.72 
+0.49 
-1.26 
+0.49 
-1.26 
+ 1.08 
-1.16 
+ 1.24 
-0.95 
+6.6 
-5.3 



4.72 

4.77 
4.43 
4.15 
5.22 
7.17 
126 
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Table 4. Sensitivity of the neutrino fluxes, the central density, pc, and the surface helium abundance, Vconv, to the uncertainties 
in the input physics.* 









CI 


Ga 


«B 




Pc 


Y 

^ conv 








7.17 


126 


4.77 




156 


0.246 


Input physics 






(SNU) 


(SNU) 


(lO^cm-^g- 




(g cm^"^) 




PP 


4.oo(i+i;-^^^)xio-^^ 


(keV b) 


— U.boicJ 
+0.372 


— 3.1 
+ 1.8 


— 0.4/7 
+0.282 




— 3.0 
+ 1.6 


+0.00U6 
-0.0003 


pep 


±1% 




±0.002 


0.0 


0.000 




0.0 


0.0000 




5.4±0.32 


(MeV b) 


TO. 131 


TO.8 


TO.094 




±0.3 


TO. 0001 




0.53±0.05 


(keV b) 


±0.408 


±2.6 


±0.293 




Tl.O 


±0.0003 


^Be+e 


±2% (la) 




TO. 106 


TO. 2 


TO.093 




0.0 


0.0000 


^Be+p 


191^(1^)^8(3^) 


(eV b) 


+ 1.057 
-0.456 


+2.2 
-1.0 


+0.927 
-0.400 




0.0 


0.0000 


i^C+p 


1.34±0.21 


(keV b) 


0.000 


0.0 


0.000 




0.0 


0.0000 


13C+P 


7.6±1 


(kcV b) 


0.000 


0.0 


0.000 




0.0 


0.0000 


"N+p 


3.5 l?:^(la)li;°(3a) 


(kcV b) 


+0.023 
-0.069 


+0.5 
-1.6 


-0.007 
+0.021 




-0.2 
+0.6 


+0.0001 
-0.0002 


^^O+p 


9.4±1.7 


(keV b) 


±0.001 


0.0 


0.000 




0.0 


0.0000 


{Zj X)surf 


0.0245±0.0006 




±0.120 


±0.5 


±0.094 




TO.2 


±0.0026 




3.844(l±0.004)xl033 


(erg s-i) 


±0.203 


±1.3 


±0.153 




±0.5 


±0.0005 




Sound-speed profile 




±0.030 


±0.1 


±0.022 




±0.3 


±0.0001 


fcomr/ Rq 


0.713±0.001 




TO.031 


TO. 2 


TO. 023 




0.0 


±0.0031 


Neutrino cross-section 




±0.173 


+5.1 
-2.4 










Opacity 


±5% 




(±)0.083 


(±)0.4 


(±)0.062 




(±)0.2 


(±)0.0075 


EOS 


ideal (r = 1) 




-0.007 


-0.3 


+0.006 




-0.6 


+0.012 


Total 


+ 0-2 




+ 1.24 

-0.95 


+6.6 


+ 1.04 
-0.72 




+2.0 
-3.2 


±0.004 


(except for opacity and EOS) 




(SNU) 


(s'nu) 


(lO'^cm-^s- 




(g cm"'^) 





* The first ten entries are due to the nuclear cross-section factors. References for each uncertainty are the same as 
given in table 1 except for (Z/X)surf- The uncertainty of {Z/X)sm:i is taken from Basu (1998). 



Table 5. Comparison of experiments and some solar models. 





Experiments 


This work 


TS98a 


TS98b 


BP98 


Cl(SNU) 


2.56 ± 0.23 * 


7 17 +1.24 

-0.95 


7.8 - 


' 10.6 


7.7 - 


- 8.8 


7 7 +^-^ 

'■' -1.0 


Ga(SNU) 


67.2 t 
77.5 tl'l * 


126 11% 


129 r 


- 142 


132 - 


- 138 


129 1% 


8B(106 cni-2 s-i) 


2.80 ± 0.38 § 

2 49 +0-12 II 

^•^^ -0.09 


A 77 

^- ' ' -0.72 


5.48 r 


- 7.67 


5.3 - 


- 6.1 


c 1 c; +0.98 
^■^"J -0.72 


Pc(g cm-3) 


153.9 ± 1.1 « 


156 t^-.i; 


156 r 


- 171 


153 - 


- 155 


-152.2 




0.248 ± 0.001 » 


0.246 ± 0.004 


0.227 r 


- 0.236 


0.246 - 


- 0.247 


-0.243 



* Homestake : Cleveland et al. 1998 

t SAGE : Abdurashitov et al. 1999 

t GALLEX : Hampel et al. 1999 

§ Kamiokande : Fukuda et al. 1996 

II Super-Kamiokande : Fukuda et al. 1998 

« Basu 1998 
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Table 6. Properties of the low-Z core model in the case of Zc = 0.0001.* 
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rf/i?0 0.000 0.050 0.069 0.083 0.10 0.12 0.13 0.15 

Mcorc/Mo 0.000 0.013 0.030 0.050 0.080 0.117 0.150 0.202 

Tc(10' K) LSe L53 L51 09 OS 06 06 05~ 

Pc (10" dyn cm-2 ) . . 2.38 2.40 2.41 2.43 2.44 2.46 2.48 2.50 

Pc (g cm-3 ) 156 157 158 159 160 162 163 164 

Xc 0.338 0.356 0.369 0.378 0.387 0.395 0.399 0.403 

Yc 0-644 0.644 0.631 0.622 0.613 0.605 0.601 0.597 

^conv (10'' K) 2^20 2l9 2l9 2^20 2^20 2l9 2l9 2l9~ 

-Pconv (10" dyn cm-2 ) 5.71 5.70 5.70 5.69 5.68 5.68 5.67 5.66 

Pconv (g cm-3 ) 0.190 0.190 0.190 0.189 0.189 0.189 0.189 0.189 

Xconv 0.736 0.737 0.737 0.736 0.737 0.737 0.737 0.738 

Fconv 0.246 0.245 0.245 0.246 0.245 0.245 0.245 0.244 

Zconv 0.0180 0.0180 0.0180 0.0181 0.0181 0.0181 0.0181 0.0181 

Cl(SNU) 7.17 6.50 5.57 4.76 3.84 3.15 2.78 2.45 

Ga(SNU) 126 122 117 114 109 105 103 101 

^6(10^ cm-^ s~M .... 4.77 4.33 3.63 3.00 2.30 1.79 1.52 1.29 



* The first two rows indicate the radius fraction and the mass fraction of the low-Z core, respectively. 



Table 7. Sensitivity of the neutrino fluxes and the central density, pc, of the low-Z core model with Zc = 0.0001 and rf/i?Q = 0.12 
{Mcoic/Mq = 0.117) to the uncertainties in the input physics.* 









CI 


Ga 


«B 


Pc 








3.15 


105 


1.79 


162 


Input physics 






(SNU) 


(SNU) 


(lOW-^s-i) 


(g cmT^) 


PP 


4.00(l+ro^3)xl0-- 


(keV b) 


-0.249 
+0.141 


-1.5 
+0.8 


-0.179 
+0.103 


-3.4 
+ 1.8 


pep 


±1% 




±0.002 


0.0 


0.000 


0.0 


^He^He 


5.4±0.32 


(MeV b) 


TO.066 


TO.6 


TO.040 


±0.3 


^He^He 


0.53±0.05 


(keV b) 


±0.207 


±1.9 


±0.125 


+0.8 


^Be+e 


±2% (la) 




TO.040 


TO.1 


TO.035 


0.0 


^Be+p 


19^4(la)^8(3a) 


(eV b) 


+0.369 
-0.171 


+0.8 
-0.4 


+0.347 
-0.150 


0.0 


i^C+p 


1.34±0.21 


(keV b) 


0.000 


0.0 


0.000 


0.0 


"C+p 


7.6±1 


(keV b) 


0.000 


0.0 


0.000 


0.0 


14N+P 


3.5 t?1(la)ii:°(3cT) 


(keV b) 


+0.000 
-0.001 


0.0 


-0.000 
+0.001 


0.0 


i^O+p 


9.4±1.7 


(keV b) 


0.000 


0.0 


0.000 


0.0 


(Z/X)surf 


0.0245±0.0006 




±0.024 


±0.2 


±0.017 


+0.2 




3.844(l±0.004)xl033 


(erg s"^) 


±0.082 


±0.8 


±0.059 


±0.6 




Sound-speed profile 




±0.011 


±0.1 


±0.008 


±0.3 


Tcomr/ Rq 


0.713±0.001 




TO.005 


TO.1 


TO.003 


±0.1 


Neutrino cross-section 




±0.066 


+3.1 
-1.7 






Opacity 


±5% 




(±)0.025 


(±)0.2 


(±)0.018 


(±)0.1 


EOS 


ideal (r = f) 




+0.001 


-0.2 


+0.006 


-1.1 


Total 


Vcr^ -hcr2 _^ 




+0.49 
-0.39 


+4.0 


+0.39 
-0.28 


+2.1 
-3.6 


(except for opacity and EOS) 




(SNU) 


(s'nu) 


(10*^cm-2s-i) 


(g cm^^) 



* The first ten entries are due to the nuclear cross-section factors. The references for each uncertainty are the same 
as shown in table 4. 
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z 







(a) 



Z 



■^urf 

Zc 



(C) 



Z 







(b) 



(d) 



+^r/R, 



Fig. 1. Comparison of Z distributions, (a) Homogeneous case (section 3, TS98a). (b) Low-Z core (section 4). (c) Evolutionary 
model with gravitational settling (e.g., BP98). (d) Rectilinear profile (TS98b). 



This work 



TS98b 




TS98; 



23456789 10 11 
V capture rate for CI experiment (SNU) 



60 80 100 120 140 

y capture rate for Ga experiment (SNU) 



This work 
TS98b 
TS98a 
BP98 

Super K 

(1998) 

Kamiokande 
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TS98a 



BP98 



(d) 
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'*B V flux(10'cm"V') 



0.23 



0.24 



0.25 



Fig. 2. Comparison of experiments and some solar models. We should regard the errors of TS98a and TS98b as rough estimations, 
(a) Neutrino capture rate for the chlorine detector at Homestakc. (b) Neutrino capture rate for the gallium detectors (GALLEX and 
SAGE), (c) *B-neutrino flux for the Kamiokande and the Super-Kamiokande. (d) Helium abundance at the base of the convection 
zone and the surface yconv. For a comparison, the scismically inverted abundance, Vconv (Basu 1998), is also shown. 




Fig. 3. Relative differences in the squared sound-speed (a) and in the density (b) between the solar models and the seismically 
determined profiles (Basu 1998). For a comparison, the profiles of an evolutionary solar model (BP98) are also shown. 



0.015 
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-a 
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(a) 
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■ +1 (J (Basu 1998) 




- This work with +1(T en'or bar (t") 




- BP98 



0.1 0.2 0.3 0.4 0.5 



0.6 0.7 




±1 a (Basu 1998) 
This work with +1 cr error bar(5| ,) 
BP98 



0.1 



0.2 0.3 0.4 0.5 0.6 0.7 
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±1 (J (Basu 1998) 
This work with ±\o en'or bar(iS34) 
BP98 
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Fig. 4. Uncertainties in the density profile caused only by errors in the sound-speed profile (a), by uncertainties in the nuclear 
cross-section 511(0) (b), and 534(0) (c). 





Fig. 5. Dependence of the neutrino fluxes and tlie central density upon Zc and the core size, (a) Neutrino capture rate for the 
chlorine detector at Homestake. (b) Neutrino capture rate for the gaUium detectors (GALLEX and SAGE), (c) *B-neutrino flux for 
the Kamiokande and the Super-Kamiokande. (d) Central density. For a comparison, the seismically inverted density, Pobs (Basu 
1998), is also shown with the observational error. 



